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ABSTRACT 


For  a long-range  radar  system  the  target  range  and  bearing 9 as  well  as 
the  signal  amplitude  and  phase.,  are  altered  to  some  degree  by  the  presence  of 
the  earth’s  atmosphere  and  ionosphere.  The  irregular  nature  of  the  ionosphere 
controls  the  extent  to  which  these  errors  can  be  predicted  and  compensated  for. 
The  limits  on  radar  metric  accuracy  imposed  by  ionospheric  propagation  effects 
are  evaluated  with  emphasis  on  an  L-band  radar  situated  near  the  auroral  zone. 
Many  of  the  results  are  derived  from  previous  propagation  studies  conducted  at 
the  Millstone  Hill  radar  facility  by  applying  appropriate  frequency  scaling 
laws.  Estimates  are  made  of  elevation  and  range  errors  associated  with  the 
ambient  ionosphere  and  with  large-scale  ionospheric  structures  (such  as  tne 
mid-latitude  electron  density  trough  and  traveling  ionospheric  disturbances). 

The  amplitude  and  angle -of -arrival  fluctuations  produced  by  ionospheric  scintil- 
lation are  considered  and  models  constructed  which  are  applicable  for  normal 
ionospheric  conditions  and  also  for  periods  of  severe  geomagnetic  disturbance. 
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I . INTRODUCTION 

This  report  reviews  the  limitations  which  the  ionosphere  imposes  on 
the  accuracy  of  radar  measurements , with  emphasis  on  an  L-band  radar  (1000  MHz) 
operating  in  or  near  the  auroral  zone.  The  effect  of  the  troposphere  on  the 
apparent  elevation  and  range  of  a target  at  low  elevation  angles  is  well 
known  and  standard  corrections  are  applied  to  allow  for  it.  Ionospheric 
effects  on  radar  propagation  are  often  ignored  because  the  necessary  correc- 
tions are  not  as  readily  applied  as  for  the  troposphere.  Also,  a rapid  decrease 
in  the  importance  of  ionospheric  errors  may  be  effected  by  increasing  the 
operating  frequency  of  the  radar.  In  some  applications  however,  the  ultimate 
accuracy  of  a radar  is  required  and  all  bias  errors  must  be  evaluated  if 
possible . 

The  ionosphere  is  defined  as  that  part  of  the  upper  atmosphere  which 
is  sufficiently  ionized  to  influence  the  propagation  of  radiowaves,  and  as 
such  extends  from  an  altitude  of  50  km  to  several  earth  radii.  While  many 
properties  of  this  region  are  still  not  understood,  the  propagation  of  radio- 
waves is  governed  solely  by  the  distribution  of  electrons  along  the  propagation 
path.  Examples  of  the  vertical  distribution  of  ionization  are  given  in 
Figure  1 for  typical  dayisime  and  nighttime  conditions  at  sunspot  maximum. 

In  practice,  the  propagation  effects  produced  by  the  ionosphere  are 
difficult  to  predict  because  of  the  large  daily,  seasonal,  sunspot-cycle, 
geomagnetic  activity  and  latitudinal  variations  in  the  ambient  electron  density. 
An  additional  source  of  propagation  errors  arises  from  the  presence  of  inhomo- 
geneities in  the  ionosphere.  These  inhomogeneities  have  a variety  of  scales 
ranging  from  global  down  to  cm  sizes. 

Many  of  the  results  reported  here  were  based  on  propagation  measure- 
ments made  at  Millstone  Hill.  These  measurements  were  made  at  UHF  and  VHF 
and  have  been  scaled  to  L-band  for  the  present  study.  However,  the  frequency 
scaling  laws  are  given,  and  with  some  restrictions,  the  results  could  readily  be 
scaled  to  other  radar  frequencies.  A brief  description  of  the  Millstone  Hill 
propagation  study  is  given  in  Section  II. 
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Fig.  1 An  idealized  electron  density  distribution  in  the  earth's  iono- 

sphere. The  curves  drawn  indicate  the  densities  to  be  expected  at  sunspot 
maximum  in  temperate  latitudes. 
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For  the  purpose  of  this  report 3 the  ionosphere  is  considered  as  a 
superposition  of  two  components;  an  ambient  unperturbed  part  and  an  irregular 
part  comprising  inhomogeneities  of  various  sizes.  The  ambient  ionosphere  is 
assumed  to  be  spherically  symmetric  in  the  vicinity  of  the  radar  and  the 
resulting  bias  errors  in  elevation  and  range  are  discussed  in  Sections  III 
and  IV.  The  propagation  effects  produced  by  various  types  of  ionospheric 
inhomogeneities  are  then  discussed  commencing  with  the  largest  scale 
phenomenon.  Section  V treats  radar  errors  at  high  latitude  which  are  associated 
with  the  mid-latitude  trough  in  electron  density  and  Section  VI  examines  the 
effects  produced  by  traveling  ionospheric  disturbances.  Sunrise  and  sunset 
effects  are  usually  negligible  and  are  ignored  in  the  present  study.  The 
effects  of  scintillations  are  treated  in  Section  VII  and  a summary  of  the 
magnitude  of  the  various  propagation  errors  at  L-band  is  given  in  the  final 
section . 
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II.  THE  MILLSTONE  HILL  PROPAGATION  STUDY 

A.  General 

Many  of  the  results  presented  in  this  report  are  based  on  a 
propagation  study  conducted  from  1969  through  1973  at  the  Millstone  Hill 
Radar  Facility  (42.6°N,  71*5°W)  jointly  by  the  staff  of  Lincoln  Laboratory 
and  Bell  Telephone  Laboratories^.  The  study  was  supported  by  the  U.S.  Army 
Advanced  Ballistic  Missile  Defense  Agency  (ABMDA)  and  the  U.S.  Army  SAFEGUARD 
System  Command. 

During  the  period  January  1971  to  March  1973 5 five  satellites  of 
the  Navy  Navigation  System  were  tracked  with  a passive  multi-lobe  UHF  beacon 
tracker  on  the  84-foot  fully- steerable  antenna.  A frequency  selective  sub- 
reflector was  mounted  in  the  antenna  and  this  allowed  simultaneous  active 
tracking  of  the  satellites  with  the  L-band  radar.  The  satellite  beacon  signal 
at  VHF  (150  MHz)  was  also  measured  with  a yagi  antenna  mounted  on  the  edge  of 
the  dish  and  combined  with  the  phase-locked  UHF  (400  MHz)  signal  to  determine 
differential-doppler . In  addition  to  the  satellite  tracking  program, regular 
incoherent  scatter  soundings  were  made  to  study  the  ionosphere. 

A number  of  scintillation  and  refraction  studies  were  undertaken 
using  the  tracking  data  and  a brief  outline  of  the  various  investigations 
which  are  relevant  to  the  present  work  is  given  here. 

B.  Scintillation  Study 

The  Navy  Navigation  satellites  were  in  almost  circular  polar 
orbits  at  heights  near  1100  km,  and  were  seen  rising  or  setting  through  the 
auroral  zone  to  the  north  where  scintillation  effects  are  often  observed.  A 
total  of  2376  satellite  tracks  were  made  during  the  measurement  period,  repre- 
senting almost  600  hours  of  observations.  The  raw  data  from  each  track  were 
recorded  on  magnetic  tape  as  digital  samples  of  various  parameters  every 
l/l5  sec. 

The  UHF  beacon  tracker  determined  the  signal  amplitude  and  the 
traverse  (traverse  = azimuth  x cosine  of  elevation)  and  elevation  angles.  A 
parameter  which  measured  the  amount  of  fluctuation  produced  by  ionospheric 
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scintillation  over  a three -second  interval  was  calculated  for  each  of  these 
three  variables. 

The  parameter  used  to  characterize  the  UHF  amplitude  fluctuations 
was  the  scintillation  index  S,  defined  as  the  rms  power  about  the  mean  divided 
by  the  mean  power 
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where  A is  the  signal  amplitude.  This  definition  of  S corresponds  to  the  Sj 

2 1 
index  of  Briggs  and  Parkin  . 

The  angle -of -arrival  fluctuations  were  measured  by  the  rms 
traverse  angle  <7  and  the  rms  elevation  angle  cr^  . In  each  case  the  rms 
values  were  calculated  after  removing  a linear  trend  by  fitting  a straight 
line  to  tie  45  sample  points.  The  trend  correction  allowed  for  relatively 
slow  changes,  which  were  present  in  addition  to  the  more  rapid  fluctuations 
produced  by  ionospheric  scintillations . 

Figure  2 is  an  example  of  the  plots  produced  by  the  post-pro- 
cessing  program,  which  calculates  three-second  values  of  track  parameters. 

The  geographic  and  invariant  latitude  scales  are  subionospheric  latitudes 
assuming  an  ionospheric  penetration  point  at  a height  of  300  km.  Millstone 
Hill  is  at  an  invariant  latitude  of  56°N. 

In  the  condensed  form  (three-second  averages)  all  the  results 
for  the  2376  satellite  tracks  were  stored  on  two  magnetic  tapes,  so  greatly 
facilitating  a statistical  analysis  of  the  scintillation  data.  Values  of  the 
geomagnetic  disturbance  index  Kp  and  local  vertical  incidence  ionosonde 
measurements  of  the  F- region  critical  frequency  fQF2  for  each  track  were 
included  on  the  condensed  data  tapes. 

A small  level  of  signal  amplitude  fluctuation  was  evident  durin 
all  satellite  tracks  which  is  not  caused  by  ionospheric  scintillation.  This 
produces  a "baseline"  level  for  the  measured  scintillation  parameters,  which 
is  evident  in  Figure  2.  The  baseline  level  is  elevation  angle  dependent  and 
is  typically  0.05  for  S.  The  two  major  sources  of  amplitude  fluctuations 
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Fig.  2 Typical  results  for  a satellite  track  obtained  for  the  Millstone  Hill 
propagation  study.  Starting  from  the  top  are  (a)  UHF-VHF  differential  phase 
<J>,  (b)  differential-doppler  d$/dt,  (c)  UHF  elevation  deviations  (3-sec  vs  63- 
sec  running  means),  (d)  UHF  traverse  deviations  (3-sec  vs  63-sec  running  means), 
(e)  fog^o  Peak  el*ectron  density  (f)  differential- phase  standard  deviation, 

(g)  UHF  elevation  standard  deviation,  (h)  UHF  traverse  standard  deviation, 
and  (i)  UHF  scintillation  index. 
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causing  the  baseline  level  were  receiver  noise  and  modulation  on  the  satellite 
beacon  signals,  although  tropospheric  scintillation  may  also  have  been  present 
at  low  elevation  angles  on  some  occasions.  A model  of  the  baseline  variation 
was  constructed  and  subtracted  from  the  measured  scintillation  levels  to  give 
the  scintillation  index  produced  by  ionospheric  irregularities  alone. 

C . Differential-Dopr ler  Experiment 

The  UHF  and  VHP  beacon  signals  were  combined  to  determine  the 
differential-doppler  d<5/dt. 


d$ 

dt 


|~VHF  phase  angle  - 3/8  UHF  phase  anglej 


(2) 


Employing  the  differential-doppler  effectively  eliminates  doppler  cnanges 
produced  by  the  changing  path  length  as  the  satellite  moves,  thereby,  leaving 
only  the  changes  produced  by  the  ionosphere. 

The  differential-doppler  d$/dt  was  essentially  proportional  to 

the  rate  of  change  of  electron  content  along  the  propagation  path  from  the 

satellite  to  antenna.  A method  of  integrating  this  measurement  to  determine 

the  variation  of  the  differential  phase  <£  and  the  peak  electron  density 

along  the  satellite  track  is  described  in  detail  by  Evans  and  Holt^.  Electron 

density  profiles  of  the  Chapman  type  were  assumed  at  all  latitudes  and 

normalization  was  achieved  using  the  local  determination  of  N from  ionosonde 

m 

measurements.  Figure  2 includes  a plot  of  log ^ derived  in  this  manner. 

The  differential-doppler  records  were  used  directly  to  give 
information  on  the  occurrence  of  long -wavelength  traveling  ionospheric 
disturbances  which  can  produce  refraction  and  range  errors.  The  derived 
parameter  permitted  a study  of  the  location  of  the  mid-latitude  electron 
density  trough  which  also  can  produce  anomalous  propagation  effects.  The 
differential  phase  $ was  especially  useful  in  this  propagation  study  as 
it  is  proportional  to  the  ionospheric  range  error. 
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D.  Refraction  Measurements 


The  L-band  radar  tracked  the  Navy  Navigation  satellites  and 
provided  a position  estimate  which  was  independent  of  that  provided  by  the 
400  MHz  beacon  tracker*  Ionospheric  refraction  at  L-band  is  negligible  in 
comparison  to  that  at  400  MHz  and  the  two  sets  of  measurements  could  be 
used  to  direct^  determine  the  gross  ionospheric  refraction  at  400  MHz. 

To  this  end*  tracking  exercises  were  conducted  over  intervals  of 
four -five  days  to  gather  sufficient  L-band  data  to  solve  for  the  orbit  of  each 
satellite  as  well  as  certain  instrumental  radar  biases.  The  true  satellite 
positions  established  in  this  manner  were  then  compared  with  the  positions 
given  by  the  beacon  tracker.  A nighttime  comparison  (when  ionospheric 
refraction  is  negligible)  allowed  further  instrumental  effects  present  in 
the  beacon  tracker  to  be  evaluated.  Finally,  the  daytime  passes  were  analyzed 
to  determine  the  refraction  produced  by  the  ionosphere  at ,400  MHz. 

A separate  investigation  of  UHF  refraction  errors  was  undertaken 
by  comparing  the  3-second  average  elevation  (and  traverse)  angle  with  the 
63-second  mean  (Figure  2).  The  procedure  eliminated  the  need  for  orbit 
determinations  and  was  suitable  for  investigations  of  the  refraction  errors 
produced  by  short -wavelength  traveling  ionospheric  disturbances. 

III.  ELEVATION  ERRORS  PRODUCED  BY  THE  AMBIENT  IONOSPHERE 

A radiowave  propagating  through  the  ionosphere  and  troposphere  under- 
goes bending  because  of  the  change  in  refractive  index  (n)  along  the  path. 

In  the  ambient  ionosphere  vertical  gradients  in  refractive  index  are  normally 
much  larger  than  any  horizontal  ones.  Consequently,  the  primary  ionospheric 
effects  may  be  estimated  by  assuming  spherical  symmetry  which  results  in  (n) 
being  a function  of  altitude  only.  Some  departures  from  the  case  of  spherical 
symmetry  are  discussed  in  later  sections. 

With  the  refractive  index  a function  of  altitude  only,  Snell ?s 
law  takes  the  form 

n r cos  E = n r cosE  (3) 
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and  the  total  tending  Acv  (Figure  3)  is  expressed  as 
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2 2 21 
n r cos  E J 
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where  Eq  and  n^  are  the  apparent  elevation  angle  and  refractive  index  at  the 
earth *s  surface , h^  is  the  height  of  the  target  and  r is  the  earth's  radius. 

The  total  bending  is  expressed  as  the  sum  of  the  ionospheric  and 
tropospheric  components  9 


Acv  = Ac^ 


+ 


(5) 


In  the  troposphere  n > 1 and  dn/dh  is  negative  causing  £&  to  be  positive 
and  the  apparent  elevation  to  exceed  the  true  elevation.  The  refractive 
index  of  the  ionosphere  for  radiowave  frequencies  (f)  in  the  VI IF  range  or 
higher  is  approximately  given  by 


n = 


1 
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4rl 


e2N 

e mf2 
o 


(6) 


where  N is  the  electron  concentration,  the  permittivity  of  free  space,  e 

the  electron  charge,  m the  electron  mass  and  f the  plasma  frequency.  The 

P 

electron  density  N is  always  sna.ll  below  60  km,  reaches  a maximum  at  height 
hm  (around  250-350  km)  and  decreases  above  h (Figure  l).  For  a radar  target 
above  a plane  ionosphere  there  would  be  no  ionosphere  bending  since  the  ray 
path  is  from  a region  where  n = 1 above  the  ionosphere,  to  a region  between 
the  ionosphere  and  troposphere  where  n is  again  unity.  However,  for  a 
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Fig.  3 Complete  ray-path  through  the  troposphere  and  ionosphere. 
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curved  ionosphere  the  upward  bending  above  h (where  dn/dh  is  positive)  is  more 
than  compensated  by  the  downward  bending  below  h (where  dn/dh  is  negative). 
Thus,  the  refraction  produced  by  both  the  ionosphere  and  troposphere  is  in  the 
same  sense  and  increases  the  apparent  elevation  angle  over  the  true  elevation 
angle. 

In  contrast  to  tropospheric  refraction,  there  is  not  a rapid  increase 
in  ionospheric  refraction  with  decreasing  elevation  angle.  This  arises  because 
the  height  of  the  ionosphere  is  non-negligible  in  comparison  with  the  earth1 s 
radius,  and  there  is  no  way  for  a ray  from  the  earth  to  arrive  at  the  iono- 
sphere at  grazing  incidence.  As  a result,  ionospheric  refraction  increases 
slowly  as  the  apparent  elevation  angle  approaches  zero  and  becomes  relatively 
constant  over  the  last  few  degrees. 

The  total  bending  Au  and  the  elevation  error  AE  (Figure  3)  may  be 
calculated  once  the  refractive  index  is  specified  as  a function  of  height. 
However,  it  is  not  generally  possible  to  derive  analytical  expressions  for  a 
realistic  ionosphere,  and  a ray  tracing  approach  must  be  used  to  determine 
Au  and  AE.  Accurate  calculations  of  ionospheric  refraction  must  include  the 
tropospheric  component  of  refraction,  as  this  alters  the  angle  of  incidence 
of  the  ray  at  the  ionosphere.  The  ionospheric  component  of  elevation  angle 
error  AE^.  may  be  obtained  from 


AEi  aetot 


AE„ 


(7) 


where  AE  results  from  the  combined  ionosphere  and  troposphere  and  AEm  is 
TOT  T 

the  result  for  the  troposphere  alone. 

The  ionospheric  component  of  elevation  error  AE  has  been  investigated 
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at  BTL  using  ray  tracing  techniques  and  electron  density  profiles  of  the 

Chapman  form 

N(h)  = Nm  exp  \ (1  - z - e“Z)  (8) 
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with 


z 


andh_  is  the  scale  height  of  the  ionizable  neutral  gas  and  N is  the 
s m 

peak  electron  density  at  height  h^.  The  plasma  frequency  at  the  peak  f 

(i.e.,  the  critical  frequency  f F2)  is  related  to  N by 
' o m 

f = 8.98  N ^ (9) 

m m 

o 

where  f is  in  the  Hz  and  N in  electrons/nr . We  have  found  that  the  BTL  ray 
m m 

tracing  results  for  AE^  can  be  fitted  reasonably  well  by  an  empirical  formula 
of  the  form 

4Et  = 31.2  X 103  / 'fnV  C°sEo  [~  E02/(21-4  + 57 -0  SXP  <-V393)f]mdes 

\ f / i + exp  [(hm  - ^)/.75  hs  ] 

(10) 

where  the  apparent  elevation  angle  Eq  in  in  degrees  and  the  target  height  h^ 

is  in  km.  Equation  10  is  applicable  for  typical  ionospheric  values  of  h 

(250-350  km)  and  h (40-60  km), 
s 

Some  of  the  BTL  ray  tracing  results  for  Chapman  profiles  are 

shown  in  Figures  4 and  5 us  functions  of  target  height  and  apparent  elevation 

angle.  Also  shown  in  the  figures  is  the  empirical  fit  given  by  equation  10. 

The  range  of  elevation  angles  available  from  the  ray  tracing  program  is  quite 

limited  (Figure  5)  and  the  empirical  fit  could  be  inadequate  for  E^  > 20°. 

Ray  tracing  studies  were  also  conducted  by  BTL  using  actual 

electron  density  profiles  as  measured  by  the  Millstone  Hill  incoherent 

scatter  radar , rather  than  profiles  having  an  assumed  Chapman  form.  The 

calculations  used  h “ 1000  km  and  f = 400  MHz.  Results  were  scaled  to 

1 -2 

f = 1000  MHz  using  a f scaling  law,  and  the  diurnal  variation  of  AE^.  is 
shown  in  Figure  6 for  observations  made  on  23-24  Feb.  1970.  On  these  days 

f attained  a maximum  daytime  value  of  12  MHz  and  an  elevation  error  of  about 

m 

4.5  mdeg  would  be  expected  at  low  elevation  angles  for  f = 1000  MHz.  The 
predictions  from  the  empirical  formula  using  the  measured  diurnal  variation 
of  are  compared  with  the  ray  tracing  results  in  Figure  7?  for  selected 
elevation  angles  of  2°,  10°  and  30°.  The  agreement  is  reasonable  for  the 
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Fig.  4 Ionospheric  elevation  errors  as  a function  of  target  height  obtained 
from  ray-tracing  through  Chapman  layers  are  compared  with  the  empirical  model 
given  in  the  text. 
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Fig.  5 Ionospheric  elevation  errors  as  a function  of  apparent  elevation 
obtained  from  ray-tracing  through  Chapman  layers  are  compared  with  the  empiri- 
cal model. 
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Fig,  6 Ionospheric  elevation  error  at  1000  MHz  computer  by  ray- tracing 
through  profiles  such  as  those  shown  in  Figure  1 for  a 24-hour  period  in 
February  1970. 
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Fig.  7 Comparison  of  elevation  errors  given  in  Figure  6 with  the  predictions 
from  the  empirical  model. 
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smaller  values  of  Eq,  but  not  so  good  for  = 30°  under  daytime  conditions. 

A further  test  of  the  empirical  model  (equation  10 ) was  made  by 
comparing  it  with  elevation  error  measurements  made  at  Millstone  Hill  (Sec- 
tion II).  The  observing  frequency  was  400  MHz,  and  the  Millstone  Hill  results 
in  Figure  8 represent  the  average  for  seven  satellite  tracking  periods 
between  August  72  and  March  73^  The  satellite  height  was  1000  km  and  the 
mean  value  of  f was  7*3  MHz  (daytime  conditions).  The  predictions  from  the 
empirical  model  provided  an  excellent  match  to  the  measurements  (Figure  8). 

Based  on  equation  10  the  elevation  error  produced  by  the  iono- 
sphere can  be  estimated  for  any  operating  frequency,  elevation  angle  and 

target  height,  once  the  ionospheric  parameters  h , h and  f are  known. 

ms  m 

Elevation  errors  are  largest  in  daytime  when  typical  values  of  h^  and  b^  are 
300  km  and  SO  km,  respectively.  The  most  important  parameter  in  determining  the 
elevation  error  is  fm,  which  undergoes  large  diurnal,  seasonal  and  sunspot 
cycle  variations.  The  maximum  values  of  fm  at  mid  and  high  latitudes  occur 
during  winter  daytime  at  sunspot  maximum,  and  only  rarely  will  f exceed 
15  MHz  under  these  conditions.  Thus,  for  an  operating  frequency  of  1000  MHz 
and  a target  height  of  1000  km  a maximum  elevation  error  of  about  7 mdeg 
could  be  expected  at  low  elevation  angles. 

In  practice  f could  be  obtained  from  ionosonde  measurements 

m 

made  at  a point  beneath  the  intersection  of  the  ray-path  and  the  level  h . 
Alternatively,  predictions  of  f such  as  those  published  by  the  Institute 
for  Telecommunications  Services  (U.S.  Dept,  of  Commerce,  Boulder,  Colorado), 
can  be  used.  This  last  approach  is  more  convenient  but  less  accurate  since 
the  predictions  can  be  in  error,  especially  during  magnetically  disturbed 
periods . 

IV.  RANGE  ERRORS  IN  THE  AMBIENT  IONOSPHERE 

In  this  section  we  examine  the  magnitude  of  the  range  errors 
introduced  by  the  unperturbed  (spherically  symmetric)  ionosphere.  The  one- 
way transmission  time  t for  a radar  pulse  is  given  by 
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Fig.  8 Comparison  of  elevation  errors  measured  during  the  Millstone  Hill 
propagation  study  at  400  MHz,  and  predictions  from  the  empirical  model. 
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(11) 


where  n' is  the  group  refractive  index,  and  the  integral  is  evaluated  along 
the  ray-path  from  the  radar  to  the  target.  As  n'  is  greater  than  unity  in 
the  ionosphere,  t will  exceed  the  free-space  propagation  time  and  cause  an 

increase  in  the  apparent  range  R - ct  over  the  line-of-sight  range  R.  The 

a 

ionospheric  radar  range  error  is  expressed  as 


AR 
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The  group  refractive  index  may  be  calculated  from 


„ dn 

n + f I? 


(13) 


where  the  ionospheric  refractive  index  n for  frequencies  f above  VTEF  is  given 

by  equation  6.  It  follows  that  for  the  ionosphere  nn'  1. 

Neglecting  ray  curvature  in  the  ionosphere  and  troposphere 
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AR  = 


n'dr  - R 


8 tt  2 e m f 
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Ndr 


(14) 


Thus,  for  high  radar  frequencies  f,  the  ionospheric  range  error  is  inversely 
2 

proportional  to  f and  directly  proportional  to  the  electron  content  along 
the  line-of-sight  to  the  target. 
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The  Chapman  electron  density  profile  (equation  8)  may  be  directly 

integrated  for  vertical  incidence  to  give  a total  electron  content  proportion 
2 

al  to  f^  hg,  where  f is  given  by  equation  9*  in  a spherically  symmetric 
ionosphere  the  total  electron  content  at  other  elevation  angles  may  be  approx 
imated  by  multiplying  the  vertical  incidence  expression  by  sec  X,  where  X is 
the  zenith  angle  at  the  point  where  the  line-of-sight  intersects  the  F2  layer 
maximum  at  height  h . Thus,  for  a Chapman  layer  we  have  the  approximate 
relationship. 


Ndr  oc  f h sec  X 
m s 


The  angle  X is  given  by 


sin  x 


r cos  E 

o o_ 

r + h 
o m 


where  r is  the  earth* s radius, 
o 


(15) 


(16) 


Ray-tracing  methods  were  used  at  BTL  to  determine  ionospheric 
range  errors  for  various  target  heights  and  elevation  angles  using  electron 
density  profiles  of  the  Chapman  form  (equation  8).  As  the  BTL  calculations 
were  only  available  for  elevation  angles  < 20°,  they  were  supplemented  with 
results  from  GE  , which  covered  the  full  range  of  elevation  angles • (The  GE 
model  included  additional  ionization  at  E and  FI  region  heights  and  gave 
slightly  different  results  than  the  BTL  model.) 

A fit  was  made  to  the  two  combined  sets  of  ray-tracing  results 
to  provide  an  empirical  expression  of  the  form 


A R 


2 


1 


1.97  h sec  X 
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+ exp 
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meters 


(17) 
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where  h is  given  in  meters.  The  form  of  equation  17  is  based  on  the  approx- 
s 

imate  expression  for  a target  at  infinity  (equations  l4,  15)  together  with  an 
added  empirical  term  in  the  denominator  to  allow  for  the  dependence  on  target 
height  hT-  The  empirical  fit  given  by  equation  17  is  compared  with  some  of 
the  ray-tracing  data  in  Figures  9 and  10. 

4 

Ray-tracing  results  were  also  performed  at  BTL  with  electron 
density  profiles  measured  with  the  Millstone  Hill  incoherent  scatter  radar  on 
February  23  and  24 , 1970.  The  calculations  assumed  spherical  symmetry  and 
took  h^  = 1000  km  and  f = 400  MHz.  The  results,  scaled  for  f = 1000  MHz  us- 

-2 

ing  a1"  f law,  are  shown  in  Figure  11  as  a function  of  local  time  and  elevation 
angle.  The  maximum  value  of  f was  12  MHz  for  these  data  and  the  corresponding 
maximum  range  error  was  48  m. 

Predictions  from  the  empirical  formula  (equation  17)  based  on  the 

measured  diurnal  variation  of  f on  February  23  and  24,  1970,  are  compared 

m 

with  the  ray-tracing  results  in  Figure  12  for  elevation  angles  of  2°,  10°  and 
30°.  In  general,  a satisfactory  agreement  is  attained. 

The  empirical  expression  for  range  error  (17)  was  further  checked 
by  comparing  with  results  from  the  Millstone  Hill  satellite  tracking  program 
(Section  II ).  Although  the  ionospheric  range  errors  were  not  directly  measured 
in  this  study,  the  line-of -sight  electron  content  could  be  estimated  and  used 
to  predict  AR  at  f = 1000  MHz  via  equation  l4.  The  content  measurements  were 
based  on  VHF-UHF  differential- doppler  measurements  combined  with  local  values 
of  f at  Millstone  Hill  (Section  II).  Range  errors  calculated  in  this  manner 
are  shown  in  Figures  13,  l4  and  15,  as  functions  of  elevation  angle  looking 
to  the  south  (azimuths  90°  to  270°)  and  looking  to  the  north  (azimuths  270° 
to  90°),  for  the  local  values  of  f at  Millstone  Hill  in  three  ranges  (2  to 
4 MHz,  8 to  10  MHz  and  10  to  12  MHz).  Each  point  represents  an  average  over 
many  satellite  tracks;  the  bar  is  the  standard  deviation  of  the  range  errors 
in  the  particular  interval,  and  thus  gives  a rough  indication  of  the  amount 
of  variability  in  range  error  expected  for  a given  local  value  of  f . The 
range  errors  are  generally  larger  to  the  south  of  Millstone  Hill  because  of 
increased  electron  concentrations  to  the  south  on  average.  The  model  pre- 
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Fig.  9 Ionospheric  range  errors  as  a function  of  target  height  from  ray- 
tracing through  Chapman  layers  compared  with  the  empirical  expression  derived 
in  the  text. 
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ELEVATION  ANGLE  (deg) 

Fig.  10  Ionospheric  range  errors  as  a function  of  apparent  elevation  angle 
obtained  from  ray-tracing  through  Chapman  layers  are  compared  with  an  empiri- 
cal model . 
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Fig.  11  Ionospheric  range  error  at  1000  MHz  computed  by  ray-tracing  through 
profiles  such  as  that  shown  in  Figure  1 for  a 24 -hour  period  in  February  1970. 
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Comparison  of  elevation  errors  given  in  Figure  11  with  predictions 
mpirical  model. 
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Fig.  13  Comparison  of  range  error  predictions  from  the  empirical  model  and 
values  derived  from  differential-doppler  data  recorded  during  the  Millstone 
Hill  propagation  study.  Results  are  for  a range  of  local  f F2  values  (f  ) 
from  2 to  4 MHz. 
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Fig.  14  Comparison  of  range  error  predictions  from  the  empirical  model  and 
values  derived  from  differential-doppler  data  recorded  during  the  Millstone 
Hill  propagation  study.  Results  are  for  a range  of  local  f F2  values  (f  ) 
from  8 to  10  MHz. 
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Fig.  15  Comparison  of  range  error  predictions  from  the  empirical  model  and 
values  derived  from  differential-doppler  data  recorded  during  the  Millstone 
Hill  propagation  study.  Results  are  for  a range  of  local  f F2  values  (f  ) 
from  10  to  12  MHz.  0 m 
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dictions  are  seen  to  provide  a good  fit  to  the  data  in  Figures  13 , 14  and  15 • 
As  the  model  is  based  on  spherical  symmetry,  the  north-south  variation  clearly 
cannot  be  matched. 

Once  the  values  of  the  ionospheric  parameters  f , h and  h are 

supplied, equation  (17)  provides  a prediction  of  ionospheric  range  error  at 

any  particular  operating  frequency,  target  height  and  elevation  angle. 

Range  errors  are  of  most  concern  during  the  daytime  when  f reaches  a peak. 

If  more  precise  information  is  lacking  typical  daytime  values  of  h^  = 300  km 

and  h = 50  km  might  be  adopted  as  these  parameters  are  not  as  important  as 
s 

f in  determining  the  range  error.  The  maximum  range  errors  would  be  encoun- 
tered at  sunspot  maximum  conditions  at  high  latitudes  in  winter,  where  f 
may  reach  15  MHz.  For  an  operating  frequency  of  1000  MHz  and  a target  height 
of  1000  km  the  maximum  radar  range  error  given  by  equation  17  is  about  75  m 
at  low  elevation  angles. 

From  equation  l4  it  is  evident  that  the  ionospheric  range  error 
is  proportional  to  the  square  of  the  radar  frequency.  Accordingly,  wideband 
FM  systems  or  other  short  pulse  systems  should  take  this  into  consideration. 

For  example,  a 10 % sweep  in  frequency  will  produce  a -20%  change  in  AR.  Thus 
for  an  FM  radar  at  1000  MHz,  a 10%  sweep  will  produce  a range  dispersion  of 
about  15  m for  the  case  when  the  radar  range  error  is  ~ 75  m (which  is  the 
maximum  ever  likely  to  be  encountered). 

V.  HIGH-LATITUDE  EFFECTS 

A.  Auroral  Clutter 

A number  of  phenomena  occur  at  high  latitudes  which  can 
influence  radar  performance.  Auroral  radar  clutter  echoes  are  observed  which 
are  associated  with  a region  known  as  the  auroral  oval  which  encircles  the 
poles  between  65°  and  80°  invariant  latitude.  These  echoes  are  of  concern 
to  a radar  looking  northward  into  the  auroral  zone  as  they  may  mask  a target 
approaching  the  radar  when  the  target  is  situated  at  the  same  range  as  the 
clutter.  Auroral  clutter  is  most  serious  when  the  radar  line-of -sight  inter- 
sects the  magnetic  field  lines  in  the  E-region  near  perpendicularity.  Because 
of  the  requirements  of  this  favored  geometry  the  clutter  is  typically  quite 
limited  in  range. 
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An  L-band  radar  study  of  auroral  returns  was  made  at  Millstone 
Hill  and  is  reported  in  detail  by  Evans1.  It  was  found  that  the  mean  height 
of  auroral  echoes  was  near  110  km  with  a mean  half -power  thickness  of  < 10  km 
and  that  the  power  fell  by  3&b  when  the  suspect  angle  with  respect  to  the 
magnetic  field  was  increased  from  0 to  .4°. 

B.  Electron  Density  Trough 

Anomalous  amounts  of  ionospheric  refraction  and  range  error 
can  be  introduced  at  high  latitudes  by  the  presence  of  large  horizontal 
(principally  latitudinal)  gradients  in  the  electron  density.  A belt  of  high 
density  is  associated  with  the  auroral  oval  where  visible  auroras  are  most 
often  seen.  Equatorwards  of  the  auroral  oval,  there  is  a narrow  latitudinal 
belt  of  low  electron  density  known  as  the  mid-latitude  electron  density 
"trough".  This  feature  is  evident  in  the  log-^Q  N plot  in  Figure  2 near  an 
invariant  latitude  of  64°N. 

The  dependence  of  the  trough  location  on  time  and  magnetic 

activity  was  established  using  the  differential-doppler  measurements  made  in 

the  Millstone  Hill  satellite  tracking  study.  For  each  satellite  track,  the 

differential-doppler  record  and  local  ionograms  were  used  to  derive  the 

variation  of  peak  electron  density  (N  ) with  invariant  latitude  as  described 

in  Section  II.  Mean  values  of  N were  calculated  for  two  hour  intervals  of 

m 

local  time  and  2 steps  of  invariant  latitude  for  two  six  month  "seasons" 
and  three  K intervals.  Sample  results  are  given  in  Figures  16  and  17  as 
contour  plots,  with  the  deduced  trough  positions  also  shown.  No  indication 
of  the  daytime  position  of  the  trough  is  evident  because  it  is  then  usually 
difficult  to  identify  such  a feature.  Mean  trough  positions  for  the  three 
K intervals  are  given  in  Figure  18.  Averages  of  the  summer  and  winter 
positions  are  shown  in  Figure  l8  as  there  was  little  seasonal  variation  evident 
for  a given  level.  In  general,  the  trough  is  seen  to  move  about  2°  south- 
ward for  unit  increase  in  K . A similar  dependence  is  given  by  Rycroft  and 
Thomas^ . 

The  trough  occurs  principally  at  night  when  the  ambient 
electron  density  is  low  and  hence  will  not  usually  affect  propagation  errors 
in  a significant  way.  However,  during  magnetically  disturbed  conditions  the 
trough  may  be  evident  during  the  daytime  and  cause  important  propagation 
anomalies . 
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Fig.  16  Contours  log^N^,  as  derived  from  differential-doppler  measurements 
made  at  Millstone  Hill  1971-1973-  Results  are  for  a interval  0Q  to  1+  in 
winter  and  the  deduced  mean  position  of  the  electron  density  trough  is  shown. 


Fie.  17  Contours  log,  11  , as  derived  from  differential-doppler  measurements 
made  at  Millstone  Hill  1971-1973.  Results  are  for  a interval  2_  to  3+  in 
summer  and  the  deduced  mean  position  of  the  electron  density  trough  is  shown. 


31 


18-2-12006-1 


Kp  00  to  1+ 

Kp  2_  to  3+ 

Kp  4_  to  5+ 


Fig.  18  Mean  location  of  tne  mid-latitude  electron  density  trough  from 
Millstone  Hill  differential-doppler  measurements . 
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In  the  Millstone  Hill  propagation  study  the  trough  was  often 
recognizable  on  single  satellite  tracks  because  of  the  perturbation  introduced 
in  the  differential-doppler  records.  Figure  19a  is  an  example  of  differential - 
doppler  data  gathered  on  successive  daytime  passes  on  one  day  when  the  trough 
could  be  seen  approaching  the  radar.  The  UHF  elevation  residuals  for  these 
passes  exhibited  corresponding  perturbations , as  shown  in  Figure  19b. 

-2 

To  estimate  the  trough-induced  elevation  errors  for  an  L-band  radar  an  f 
scaling  law  was  applied  to  the  measured  UHF  elevation  errors.  Thus,  Figure  19b 
includes  a scale  appropriate  to  an  operating  frequency  of  1000  MHz  and  it  is 
seen  that  in  this  instance  the  trough  associated  refraction  errors  were  about 
+3  mdeg  with  a decay  after  sunset.  Ray-tracing  studies  have  established  that 
these  positive  elevation  excursions  are  associated  with  rays  that  traverse 
the  lower  ionosphere  equatorwards  of  the  trough  and  emerge  above  the  layer 
peak  in  the  trough,  i.e.,  where  the  densities  are  low.  A corresponding 
negative  excursion  associated  with  the  polar  edge  of  the  trough  is  not  seen, 
apparently,  because  the  horizontal  gradients  in  density  are  smaller  at  the 
poleward  edge  and  because  the  required  geometry  is  less  readily  achieved.'1' 

Based  on  the  results  of  the  Millstone  tracking  study  at  UHF  the  largest 
trough-induced  elevation  errors  at  L-band  would  have  been  about  +9  mdeg. 

Range  errors  could  also  be  estimated  from  the  differential- 
doppler  data  because  both  the  radar  range  error  and  the  differential-phase 
angle  introduced  by  the  ionosphere  are  proportional  to  the  total  electron 
content  along  the  radar  line-of-sight  to  the  satellite.  For  a radar  frequency 
f (in  MHz)  the  ionospheric  radar  range  error  (in  meters)  is  given  by, 


AR 


8340  $ 


m 


(18) 


where  $ is  the  UHF-VHF  differential-phase  angle  in  radians . To  use  this 
relationship  to  estimate  range  errors,  it  liras  first  necessary  to  integrate  the 
values  of  d<J>/dt  obtained  along  a satellite  track  and  evaluate  the  constant  of 
integration.  The  procedure  which  was  implemented  assumed  electron  density 
profiles  of  the  Chapman  type,  and  used  local  fQF2  measurements  to  determine 
the  peak  electron  density  (Section  II).  Figure  19c  gives  the  range  errors 
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Fig.  19a  Successive  records  on  28  November  1972,  (passes  of  the  Navy  Naviga- 
tion Satellites)  exhibiting  tracking  perturbations  in  UHF-VHF  differential - 
doppler  introduced  by  the  presence  of  the  high-latitude  trough. 
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Fig.  19b  Successive  records  on  28  November  1972?  (passes  of  the  Navy  Naviga- 
tion Satellites)  exhibiting  tracking  perturbations  in  UHF  elevation  an^le 
introduced  by  the  presence  of  the  high-latitude  trough.  Dashed  line  marks 
the  location  of  the  peak  differential-doppler  excursion. 
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IONOSPHERIC  RANGE  ERROR  AT  1000  MHz  (m) 


Fig.  19c  Successive  records  on  28  November  1972 > (passes  of  the  Navy  Naviga- 
tion Satellites)  showing  the  range  error  at  1000  MHz  introduced  by  the 
presence  of  the  high-latitude  trough.  Dashed  line  marks  the  location  of  the 
peak  differential-doppler  excursion. 
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at  f = 1000  MHz  derived  from  equation  18  and  the  differential-doppler  results 
in  Figure  19a.  For  this  example  the  trough  has  caused  range  errors  about  5 m 
less  than  would  have  been  produced  by  the  ambient  ionosphere  in  the  absence 
of  the  trough.  From  the  largest  trough-induced  differential-doppler  perturba- 
tion observed  during  the  Millstone  study,  we  estimate  the  maximum  range  error 
would  be  about  -10  m. 

VI.  TID  - INDUCED  ERRORS 

Large  scale  electron  density  perturbations  in  the  F-region  of  the 
ionosphere  have  been  observed  by  a number  of  different  techniques.  These 
wavelike  electron  density  fluctuations  characteristically  possess  a horizontal 
component  of  velocity  and  are  referred  to  as  traveling  ionospheric  disturbances 
(TIDs).  TIDs  are  believed  to  result  from  internal  atmospheric  gravity 
waves  which  propagate  in  the  neutral  gas  and  produce  compression  and  rare- 
faction of  the  electron  gas.  Wave  periods  from  minutes  to  hours  and  horizontal 
wavelengths  from  tens  to  thousands  of  kilometers  have  been  observed.  The 
slowly -varying  angle  and  range  biases  produced  by  TIDs  can  degrade  a radar 
system fs  ability  to  locate  a target  and  predict  its  future  position. 

It  is  expected  that  the  approximations  of  geometrical  optics  will 
apply  for  propagation  through  TIDs  so  that  angle  and  range  errors  will  scale 
inversely  as  radar  frequency  ( f ) squared.  For  the  purpose  of  this  study, 
free  use  is  made  of  TID  results  obtained  at  UHF  by  applying  the  f scalin 
factor  to  extrapolate  to  L-band. 

The  experimental  study  conducted  at  Millstone  Hill  (Section  II ) 

demonstrated  some  of  the  propagation  effects  produced  by  TIDs.  TID-induced 
fluctuations  in  differential-doppler  and  elevation  angle  were  evident  in  many 
of  the  satellite  tracks  observed  at  Millstone  Hill,  and  an  example  is  shown 
in  Figure  20.  The  fluctuations  in  traverse  angle  were  generally  smaller  than 
in  elevation  angle  (Figure  20 ). 
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Fig.  20  Example  of  correlated  differential-doppler  and  elevation  error  per- 
turbations at  UHF  obtained  from  the  Millstone  Hill  propagation  study. 
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A statistical  study  of  the  occurrence  and  amplitude  of  TIDs  was 
based  on  the  differential-doppler  data  recorded  over  the  two  year  period  of 
the  Millstone  Hill  satellite  tracking  program.  The  rapid  motion  of  the 
satellite  ensured  that  the  time  for  the  radar  line-of-sight  to  traverse 
a TID  wavelength  was  shorter  than  the  TID  period,  so  that  the  differential- 
doppler  records  revealed  the  horizontal  spatial  structure  of  a TID,  but  gave 
no  information  on  its  period.  For  each  TID,  the  wavelength,  amplitude  and 
latitude  at  the  peak  differential-doppler  excursion  were  scaled  from  the 
tracking  records,  and  used  to  construct  the  histograms  shown  in  Figures  21, 

22,  23  and  2k. 

There  was  a marked  diurnal  variation  in  the  occurrence  of  TIDs 
with  a maximum  during  the  daytime  (Figure  21).  A partial  explanation  of  this 
effect  is  that  the  electron  density  fluctuations  produced  by  a given  gravity 
wave  are  proportional  to  the  ambient  electron  density,  which  is  smaller  at 
night . 

The  geographic  distribution  of  TIDs  (Figure  22)  shows  a strong 
peak  to  the  south  of  Millstone  Hill,  although  some  TIDs  are  observed  to  the 
north.  The  phase  fronrs  associated  with  a TID  are  tilted  in  the  direction  of 
travel  and  anomalous  propagation  effects  are  most  pronounced  when  the  direction 
of  propagation  is  aligned  along  a phase  front.  Thus,  the  peak  in  occurrence  to 
the  south  is  interpreted  as  evidence  that  most  TIDs  seen  at  Millstone  Hill  are 
produced  by  southward  propagating  gravity  waves.  This  indicates  a generation 
mechanism  to  the  north,  possibly  in  the  auroral  zone.  TID-induced  angle 
fluctuations  were  observed  mainly  in  elevation  rather  than  in  azimuth,  which 
is  further  evidence  that  the  TID  wavefronts  are  oriented  E-W  rather  than  N-S. 

Figure  23  shows  a peak  in  the  wavelength  distribution  of  TIDs 
near  250  km.  The  detailed  shape  of  the  wavelength  distribution  may  be  partly 
due  to  observational  biases,  but  the  rapid  fall  off  at  short  wavelengths  is 
apparently  a true  effect  for  differential-doppler  fluctuations  produced  by 
TIDs.  However,  TIDs  of  shorter  wavelengths  are  still  present  and  they  can 
produce  elevation  angle  fluctuatuions , as  discussed  later.  The  wavelength 
distribution  in  Figure  23  can  be  reasonably  represented  as, 
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Fig.  21  Percentage  of  tracks  exhibiting  TIDs  in  the  differential-doppler 
records  gathered  at  Millstone  Hill  between  January  1971  and  March  1973  vs  time 
of  day. 
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ig.  22  Number  of  tracks  exhibiting  TIDs  in  the  differential-doppler  records 
athered  at  Millstone  Hill  between  January  1971  and  March  1973  vs  latitude. 
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Fig.  23  Number  of  tracks  exhibiting  TIDs  in  the  differential -doppler  records 
gathered  at  Millstone  Hill  between  January  1971  and  March  1973  vs  TID  wavelength. 
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Fig.  2k  Number  of  tracks  exhibiting  TIDs  in  the  differential-doppler  records 
gathered  at  Millstone  Hill  between  January  1971  and  March  1973  vs  TID  amplitude. 
The  amplitude  has  been  converted  to  equivalent  elevation  angle  amplitude  at  UHF 
(peak-to-peak)  using  a scaling  law  of  10  mdeg/Hz. 


43 


(19) 


% of  TIDs  with  wavelength  exceeding  X (km) 

X 


= 100  exp 
£ 100  km 


-(X/350)2 


The  satellite  tracking  data  were  also  analyzed  to  determine  the 
ratio  of  the  amplitudes  of  the  TID-induced  fluctuations  in  elevation  angle 
and  differential-doppler . The  resulting  average  value  was  near  10  mdeg/Hz, 

] 

and  was  similar  to  values  predicted  by  ray -tracing  studies  using  TID  models. 
With  this  scaling  factor  the  measured  distribution  of  TID  amplitudes  from 
the  differential-doppler  records  was  converted  to  elevation  angle  and  scaled 
to  a frequency  of  1000  MHz.  Figure  24  gives  the  resulting  amplitude  distri- 
bution for  the  2CPfo  of  the  satellite  tracks  for  which  TIDs  were  observed; 

denotes  the  peak-to-peak  amplitude  of  the  TID-induced  fluctuation  in 
elevation  angle.  The  distribution  in  Figure  24  (for  a 1000  MHz  operating 
frequency)  is  approximately  given  by, 

% of  TIDs  with  peak-to-peak  amplitude  exceeding  A(mdeg)  = 100  exp 

X ^ 100  km 

(20) 

The  radar  range  error  and  the  differential-phase  angle  intro- 
duced by  the  ionosphere  are  both  proportional  to  the  total  electron  content 
along  the  radar  line-of -sight  to  the  satellite,  as  given  in  equation  l8.  An 
example  of  range  errors  calculated  from  differential-doppler  data  via  equa- 
tion 18  is  shown  in  Figure  25 , along  with  the  measured  UHF  elevation  error 
scaled  to  1000  MHz.  In  this  example,  there  are  TID-induced  range  errors 
of  4 m peak-to-peak  and  elevation  errors  of  10  mdeg  peak-to-peak,  which  are 

superimposed  on  the  values  produced  by  the  ambient  ionosphere. 

4 

Ray-tracing  studies  using  TID  models  give  the  approximate 
relationship  between  the  amplitudes  of  TID-induced  range  and  elevation  errors 
( in  m and  mdeg,  respectively)  as, 


[-  (A/2)2] 
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Fig.  25  Example  of  TID-induced  elevation  and  range  errors  (scaled  to  1000 
MHz)  from  Millstone  Hill  satellite  tracking  record. 
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100  < X < 2000  km 

where  X is  the  TID  wavelength  in  km. 

The  approximate  amplitude  distribution  of  range  errors  produced 
by  TIDs  may  be  obtained  by  combining  equations  (21)  and  (20).  Equation  (2l) 
indicates  that  range  errors  may  become  less  important  than  elevation  errors 
for  short  wavelength  TIDs. 

Predictive  schemes  may  be  used  to  correct  for  most  of  the  radar 
range  errors  and  angle  biases  produced  by  the  ambient  ionosphere  depending 
on  how  accurately  f is  estimated  (Section  III  and  IV).  However,  TID-induced 
errors  cannot  be  predicted  for  an  individual  radar  track  as,  at  best,  only 
probability  of  occurrence  can  be  specified.  One  real-time  corrective  scheme 
for  range  error  involves  simultaneous  radar  measurements  of  propagation  time 
at  two  different  frequencies.  The  possible  presence  of  TID-induced  range 
errors  places  a limit  on  how  reliably  range  error,  determined  in  this  way, 
can  be  extrapolated  to  other  points  along  the  track  or  to  other  nearby 
objects,  as  the  TID-induced  range  error  may  change  sign  for  penetration 
points  with  a horizontal  separation  of  X^-j.^/2.  The  problem  of  TID-induced 
errors  is  of  most  concern  in  the  daytime  for  a radar  situated  below  the 
auroral  zone  and  looking  southward. 

The  TID  characteristics  discussed  so  far  apply  to  wavelengths 
above  100  km,  as  TIDs  of  shorter  wavelengths  apparently  did  not  produce 
differential-doppler  signatures.  The  cutoff  at  short  wavelengths  displayed 
by  the  differential-doppler  data  was  not  present  in  TID-induced  elevation 
angle  fluctuations,  as  was  evident  from  the  UHF  beacon  tracking  results. 

In  this  study  the  UHF  traverse  and  elevation  angle  fluctuations  were  deter- 
mined by  comparing  three-sec  averages  with  63-sec  running  averages  along  the 
satellite  track.  The  smoothing  process  effectively  limited  the  range  of  TID 
wavelengths  which  could  be  studied  to  bet ween  20  and  150  km.  The  angle 


46 


fluctuations  detected  in  this  way  were  confined  mainly  to  elevation,  and  an 
example  is  shown  in  Figure  26,  scaled  to  an  operating  frequency  of  1000  MHz. 

The  elevation  angle  fluctuations  are  often  almost  sinusoidal 
in  character  (as  in  Figure  26),  but  sometimes  appear  to  be  the  mixture  of  a 
number  of  different  wavelengths.  In  order  to  characterize  the  TID  amplitude 
in  an  objective  way,  the  equivalent  peak-to-peak  amplitude  (AE^^)  was  cal- 
culated as  2/2  x rms  elevation  angle  fluctuation,  as  would  apply  for  a 
purely  sinusoidal  oscillation.  The  mean  TID  amplitude  (AEr  ()  obtained  by 
averaging  over  all  satellite  tracks  is  shown  as  a function  of  elevation  angle 
in  the  lower  part  of  Figure  27 , both  for  looking  south  and  looking  north. 

In  contrast  with  the  larger  wavelength  TIDs  observed  in  the  differential- 
doppler  study,  there  is  little  tendency  for  more  short  wavelength  TIDs  in  the 
south  than  in  the  north.  This  argues  for  a different  generation  mechanism 
for  these  TIDs  which  is  not  associated  with  the  auroral  zone. 

The  diurnal  variation  of  mean  1ID  amplitude  (scaled  to  1000  MHz) 
is  presented  in  the  upper  part  of  Figure  27-  The  maximum  mean  amplitude 
occurs  in  the  daytime,  as  with  longer  wavelength  TIDs,  and  the  explanation 
again  may  be  the  controlling  influence  by  the  diurnal  variation  in  ambient 
electron  density. 

The  distribution  of  the  peak-to-peak  amplitude  of  the  short 
wavelength  TIDs  is  sho wn  in  Figure  28.  All  satellite  tracks  were  included 
in  this  analysis,  rather  than  just  those  exhibiting  recognizable  TID-induced 
elevation  errors. 

VII.  IONOSPHERIC  SCINTILLATION 

A.  General 

Observations  by  means  of  a variety  of  techniques  have  shown  that 
small-scale  irregularities  in  the  ionization  concentration  (and  hence  refrac- 
tive index)  are  a normal  feature  of  the  ionosphere.  There  is  a marked 
variation  in  the  incidence  of  irregularity  structure  over  the  globe  with  the 
maximum  rate  of  occurrence  and  intensity  in  the  equatorial  zone  at  night  and 
in  the  polar  and  auroral  zones.  In  addition  to  latitudinal  and  diurnal 
variations,  the  occurrence  of  irregularities  is  correlated  with  season,  sunspot 
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Fig.  26  Example  of  elevation  and  traverse  angles  measured  by  the  Millstone 
Hill  UHF  satellite  beacon  tracker.  Angles  were  computed  as  a 3-sec  mean  value 
referenced  to  63-sec  mean  value. 
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tig.  27  Local  time  dependence  (upper  part)  and  elevation  angle  dependence 
(lower  part)  of  the  mean  peak-to-peak  amplitude  of  elevation  angle  fluctuations 
induced  by  short  wavelength  TIDs.  The  mean  was  obtained  by  averaging  over  all 
the  UHF  beacon  tracker  records  (Figure  26)  gathered  at  Millstone  Hill  from 
January  1971  to  March  1973* 
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Fig.  28  Distribution  of  peak-to-peak  amplitude  of  elevation  angle  fluctuations 
produced  by  short  wavelength  TIDs,  as  observed  in  the  Millstone  UHF  beacon 
tracker  records  gathered  from  January  1971  to  March  1973- 
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number  and  magnetic  activity. 

The  fluctuations  imposed  on  radio  waves  which  are  scattered  by 
the  irregularities  are  called  scintillations  and  are  present  in  amplitude, 
phase  and  angle -of -arrival  of  the  received  radio  wave.  The  irregularities 
causing  scintillation  of  signals  above  VHF  are  almost  entirely  confined  to 

the  F-region  and  are  elongated  along  the  geomagnetic  field.  Scintillation 

-2 

effects  tend  to  decrease  with  increasing  frequency,  but  the  simple  f 
scaling  law  appropriate  to  geometric  optics  is  not  necessarily  applicable  to 
scintillation. 

A study  of  scintillation  in  the  auroral  zone  was  undertaken 
using  the  UHF  (U00  MHz)  amplitude  and  angle  data  collected  on  over  2300 
satellite  beacon  tracks  conducted  for  the  Millstone  Hill  propagation  study 
(Section  II).  The  results  of  this  study  allowed  models  to  be  constructed 
giving  the  occurrence  of  scintillation  as  a function  of  invarient  latitude 
(U2C -70°N) , local  time,  geomagnetic  activity  (K  ) and  season,  although  the 
model  could  still  only  strictly  be  applied  to  the  sunspot  conditions  pre- 
vailing when  the  observations  were  made  (1971-1973)-  The  models  were 
extended  to  L-band  by  employing  frequency  scaling  laws  based  on  an  assumed 
three-dimensional  power-law  wavenumber  spectrum  of  electron  density  fluctua- 
tion with  an  exponent  of  -h . 

B . Elevation  Angle  Dependence  of  Scintillation 

The  manner  in  which  the  scintillation  index  varied  along  a 
satellite  track  was  a combination  of  a true  spatial  change  in  the  irregular- 
ity characteristics  and  effects  produced  by  the  changing  geometry  along  the 
track.  For  a uniform  layer  of  irregularities  at  a fixed  height,  the  scintil- 
lation index  increases  as  the  elevation  angle  is  lowered  because  of  the  in- 
creasing thickness  of  the  irregularity  layer  along  the  propagation  path,  and 
also  because  the  observer  moves  further  from  the  irregularities.  If  these 
purely  geometrical  effects  are  corrected  for,  then  the  true  latitudinal 
dependence  of  the  scintillation  index  and  irregularity  characteristics  may 
be  studied. 


The  elevation  scaling  law  for  the  scintillation  index  was  deter- 
mined directly  from  the  Millstone  Hill  satellite  tracking  results.  This  was 
possible  because  the  wide  variety  of  satellite  tracking  geometries , which 
were  used  in  the  study,  permitted  observations  at  the  same  invariant  latitude 
for  a range  of  different  elevation  angles.  In  each  2°  interval  of  invariant 
latitude,  the  mean  scintillation  index  was  determined  as  a function  of 
elevation  angle  and  used  to  calculate  the  elevation  correction  factor.  The 
results  are  shown  in  Figure  29  along  with  the  theoretical  curve  for  a power- 
law  irregularity  spectrum  with  a three-dimensional  spectral  index  of  minus 
four.  This  type  of  irregularity  spectrum  predicts  the  observed  elevation 

dependence  quite  well  and  explains  many  other  characteristics  of  ionospher- 

7 

ic  scintillation  observed  at  Millstone  Hill  and  elsewhere  . Based  on  the 
elevation  scaling  law  shown  in  Figure  29  all  the  measured  scintillation 
indices  were  converted  to  the  standard  case  of  a transmitter  at  an  infinite 
distance  in  the  zenith.  This  procedure  removed  any  dependence  of  the  results 
on  the  particular  latitude  at  which  the  observations  were  made  and  the  height 
of  the  satellite  beacon  which  was  tracked. 

The  ionospheric  irregularities  producing  scintillation  are 
elongated  along  the  earth Ts  magnetic  field  and  the  scintillation  index  depends 
on  the  elongation  ratio  a and  the  propagation  angle  $ (angle  between  propa- 
gation direction  and  magnetic  field).  Thus,  the  results  should  be  normalized 
to  some  reference  value  of  as  well  as  a reference  elevation  angle.  How- 
ever, the  dependence  of  S on  i()  and  a is  not  strong  provided  that  j/j  > 30°,  a 
condition  which  is  satisfied  by  most  of  the  Millstone  Hill  satellite  tracks, 
and  it  was  sufficiently  accurate  to  ignore  the  elongation  of  the  irregulari- 
ties in  normalizing  the  observations  of  scintillation  index.  It  should  be 
remembered  that  the  results  only  strictly  apply  for  observations  transverse 
to  the  field  lines  rather  then  along  the  field  lines. 

C . Dependence  on  Geomagnetic  Activity  and  Season 

A statistical  analysis  of  the  Millstone  Hill  tracking  results 
was  carried  out  after  correcting  the  measured  scintillation  indices  for 
baseline  effects  and  normalizing  to  a transmitter  at  infinity  in  the  zenith. 


52 


CORRECTION  FACTOR  FOR  SCINTILLATION  INDEX 


ELEVATION  ANGLE  (deg) 

Fig.  29  Elevation  dependence  of  scintillation  index  from  Millstone  Hill 
satellite  tracking  data. 
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The  data  were  divided  into  three  levels  of  geomagnetic  activity  (quiet,  K from 

P 

0 to  1 ; moderate.  K from  2 to  3 5 and  disturbed,  K from  4 to  5 ) and  two 
seasons  (summer  and  winter?  defined  between  equinoxes).  For  each  two  hour 
interval  of  local  time  and  2°  interval  of  invariant  latitude,  the  mean  scint- 
illation index  was  calculated  and  used  to  construct  contour  plots.  The  results 
are  given  in  Figures  30a  to  30f  for  the  different  levels  and  seasons. 

In  general,  scintillation  activity  is  confined  to  invariant 
latitudes  above  55°  and  increases  with  increasing  K^.  Some  quite  marked 
seasonal  differences  are  also  evident.  In  summer  the  mean  scintillation 
index  is  a maximum  around  local  time  20  to  02  hours,  whereas  in  winter  the 
maximum  occurs  around  14  to  20  hours, i.e.,  about  six  hours  earlier. 

D.  Scintillation  Boundary 

The  scintillation  boundary  is  defined  as  the  equatorward  edge  of 

F-region  irregularity  structure  in  the  auroral  zone  which  causes  scintillation 

8 9 

of  UHF-VHF  radio  waves" * . The  present  observations  on  occasion  show  a rapid 
increase  in  scintillation  at  high  latitudes  which  is  characteristic  of  a sharp 
boundary  (see  Figure  2 for  example),  but  more  often  the  increase  is  more  in 
the  nature  of  a gradual  onset  of  scintillation  activity.  A scintillation 
boundary  can  be  defined  from  the  present  statistical  analysis  in  an  objective 
way  by  designating  a particular  (but  arbitrary)  value  of  the  mean  scintilla- 
tion index  to  represent  a boundary  level.  Figures  31a  and  31b  show  the  scin- 
tillation boundaries  which  result  from  selecting  a mean  scintillation  index 
of  0.02  as  a boundary  value. 

The  seasonal  dependence  of  scintillation  is  clearly  evident  in 
the  position  of  the  scintillation  boundary  shown  in  Figures  31a  and  31b. 

There  is  a southward  progression  of  the  boundary  with  increasing  K , which 
is  especially  marked  in  the  midnight  sector  in  winter  and  the  early  morning 
sector  in  summer.  This  picture  of  the  scintillation  boundary  is  somewhat  de- 
pendent on  the  particular  mean  value  of  scintillation  index  selected  as  a 
boundary  value. 


54 


LU 

Q 

3 


< 

or 

5 


12  14  16  18  20  22  24  2 4 6 8 10  12 


LOCAL  TIME  (hours) 


Fig.  30a  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 

from  Millstone  Hill  satellites  tracking  observations  1971 “1973 > for  summer 

with  K from  0 to  1 . 
p o + 


LOCAL  TIME  (hours) 


Fig.  30b  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 
from  Millstone  Hill  satellite  tracking  observations  1971-1973}  for  summer  with 

K from  2 to  3 • 

P + 
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Fig-  30c  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 
from  Millstone  Hill  satellite  tracking  observations  1971-1973 , for  summer  with 

K from  4 to  5 • 

P + 


LOCAL  TIME  (hours) 

Fig.  30d  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 

from  Millstone  Hill  satellite  tracking  observations  1971-1973?  for  winter  with 

K from  0 to  1 , • 
p o + 
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Fig.  30e  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 
from  Millstone  Hill  satellite  tracking  observations  1971-1973}  for  winter  with 

K from  2 to  3 . 

P + 


LOCAL  TIME  (hours) 


Fig.  30f  Contours  of  mean  scintillation  index  at  400  MHz  (normalized  to  zenith) 
from  Millstone  Hill  satellite  tracking  observations  1971-1973}  for  winter  with 

K from  4 to  5 • 

P + 
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18-2-12007-1 


Kp  00  to  1+ 

Kp  2_  to  3+ 

Kp  4_  to  5+ 


Fig.  31a  Positions  of  the  scintillation  boundary  in  winter  for  three  K 
intervals.  The  boundary  is  defined  by  a mean  scintillation  index  of  0.?2  at 
400  MHz  and  vertical  incidence. 
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18-2-12008-1 | 


Kp  00  to  1 + 

Kp  2_  to  3+ 


Kp  4_  to  5+ 


Fig.  31b  Positions  of  the  scintillation  boundary  in  summer  for  three  Kp 

intervals.  The  boundary  is  defined  by  a mean  scintillation  index  of  0.02  at 
400  MHz  and  vertical  incidence. 
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The  position  of  the  trough  (Section  V)  is  also  indicated  on  the 
scintillation  contour  plots  in  Figures  30a  to  30f . No  strong  correlation 
between  the  occurrence  of  scintillation  and  the  trough  position  is  evident 
from  this  comparison. 

An  analysis  of  satellite  tracks  for  which  both  a sharp  scintilla- 
tion boundary  and  a trough  were  present  similary  showed  a lack  of  correlation 
in  the  positions  of  the  two  features.  Figure  2 provides  a good  example  of 
this  poor  correlation. 

E.  Model  for  Scintillation  Index 

The  scintillation  results  may  be  used  to  generate  a model  varia- 
tion of  scintillation  index  for  use  in  tracking  simulation  studies  where  the 
actual  signal  amplitude  (and  angle -of -arrival)  is  required  along  the  track. 

The  first  step  is  to  determine  the  variation  in  the  mean  scintillation  index 
with  invariant  latitude  S (A)  along  the  track  for  the  appropriate  local  time, 

K , and  season  (Figures  30&  to  30f).  The  irregularities  may  be  considered 

3? 

to  be  at  a height  of  300  km,  and  A is  the  invariant  latitude  of  the  subiono- 
spheric  point  at  the  penetration  point  of  the  irregularities.  The  variation 
of  S calculated  in  this  way  represents  the  average  to  be  expected  using  a 
large  number  of  tracks . Figure  32  gives  the  distribution  function  for  the 
ratio  of  scintillation  index  to  mean  scintillation  index  derived  from  the 
Millstone  Hill  tracking  data.  Thus,  a representative  track  for  a particular 
tracking  simulation  should  be  multiplied  by  a factor M 3 which  is  selected  on 
a random  basis  from  a distribution  having  the  form  shown  in  Figure  32. 
Allowance  is  made  for  the  particular  observing  frequency  f(in  MHz),  and 
elevation  angle  E using  results  applicable  to  the  assumed  power-law  irregular- 
ity spectrum  and  weak-scatter . Thus,  the  final  expression  for  the  model 
scintillation  index  is 

s (A)  = s (a )m  (-1 


sec  i 


(22) 
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Fig.  32  Probability  distribution  of  scintillation  index  divided  by  mean 
scintillation  index,  as  derived  from  the  Millstone  Hill  satellite  tracking 
study. 
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where 


z 


Z1Z2 


Z1  + Z2 


and.  the  parameters  z. z^,  h^  and.  i are  defined  in  Figure  33-  The  zenith  angle 
i at  the  penetration  point  of  the  irregularities  is  related  to  the  elevation 
angle  by 

i = arcsin 


r cos  E 

o 

r + hu 


(23) 


Equation  22  only  applies  for  z > 0,  and  S should  be  set  to  zero  once  the 
irregularity  layer  is  penetrated.  Physically  it  is  not  possible  to  have  S 
much  greater  than  unity,  so  if  equation  22  gives  a value  S > 1 then  we  should 
set  S = 1. 

F.  Model  for  RMS  Amplitude  Fluctuations 

The  received  amplitude  (A)  of  a radiowave  after  passing  through 

a layer  of  ionospheric  irregularities  will  fluctuate  about  the  mean  value  A. 

The  rms  fluctuation  in  relative  amplitude  R = A/A  is  specified  by  a , 

R 


a 


2 

R 


(24) 


The  scintillation  index  S used  in  the  measurements  made  at 
Millstone  Hill  was  based  on  the  rms  fluctuation  in  received  power  relative 
to  the  mean,  i .e . , 


S 


2 


(25) 


For  small  fluctuations  in  amplitude  we  may  convert  between  the  two  measures 
of  scintillation  using 


s = 2 aR  (26) 

Thus  the  model  value  of  scintillation  index  effectively  yields  a value  for 

a via  equation  26. 

K 
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G.  Probability  Distribution  of  Amplitude 

Theoretical  and  experimental  results  indicate  that  the  probability 
distribution  of  amplitude  fluctuations  produced  by  ionospheric  scintillation 
is  satisfactorily  represented  by  the  Nakagami  m-distribution, 


VB> 


2m 


m 


2 

^m-l  -mB 
B e 


r(m) 


B s 0 
m s ^ 


(27) 


where  B is  the  amplitude  normalized  by  dividing  by  the  rms  amplitude, 


B = 


R 


1 + a 


R 


(28) 


The  parameter  m is  related  to  an  by  m = (2  a ) 

K R 

H.  Allowance  for  Radar  Case 

Once  a value  of  cr  is  calculated  from  the  model,  the  actual 
K 

relative  amplitude  R may  be  selected  on  a random  basis  from  a distribution 
given  by  equation  27-  The  relative  amplitude  so  derived  applies  to  a one- 
way transmission  through  the  irregularities  and  must  be  adjusted  to  the  radar 
situation  which  involves  a two-way  transmission.  As  time  from  transmission 
to  reception  of  a radar  pulse  is  much  less  than  the  correlation  time  of  the 
amplitude  fluctuations,  the  conversion  to  the  radar  case  is  accomplished  by 
squaring  the  relative  amplitude.  Final  conversion  to  the  received  radar 
amplitude  (c)  is  realized  by  multiplying  by  the  mean  value  of  radar  amplitude 
(C),  i.e.,  the  amplitude  if  no  scintillations  were  present.  Thus 


C = R2  C (29) 

I.  Correlation  Time  of  Amplitude  Fluctuations 

The  fluctuations  in  amplitude  and  angle -of -arrival  at  a receiver 
arise  either  because  the  ionospheric  irregularities  are  changing  with  time, 
or  because  the  diffraction  pattern  is  moving  past  the  receiver  due  to  motion 
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of  the  source  or  the  irregularities . The  effects  due  to  change  in  the  iono- 
sphere are  usually  negligible  and  it  may  often  be  assumed  that  the  source  is 
moving  rapidly  enough  so  that  the  irregularity  motion  may  be  ignored.  Roughly 
speaking,  the  correlation  time  is  the  time  taken  for  the  penetration  point  at 
neight  h (Figure  33)  to  traverse  a "typical"  irregularity  scale  size.  In  gen- 
eral, a mixture  of  scale  sizes  are  present,  but  the  scale  size  most  effective 
in  producing  amplitude  scintillation  has  a radius  of  one  Fresnel  zone,  i.e., 

/X z,  where  \ is  the  observing  wavelength.  Thus,  in  a first  approximation, 
the  correlation  time  t for  amplitude  fluctuations  may  be  taken  as 


/Xz~ 


(30) 


\'/here  the  effective  penetration  point  speed  V is  given  by 


V = ZJL  (zx  + z2)(e2  + T2)"  (31) 


E and  T are  the  rates  of  change  of  elevation  and  traverse  angle  (traverse  - 
azimuth  x cosine  of  elevation),  respectively,  expressed  in  radians/sec. 

In  some  circumstances  the  target  may  be  moving  such  that  is 
small  and  the  correlation  times  are  determined  by  the  motion  of  the  irregular- 
ities. This  is  allowed  for  by  adopting  a minimum  value  of  100  m/sec  for  V , 

1? 

corresponding  to  a most  probable  irregularity  speed. 

J . Model  for  Angle  Fluctuation 

The  UHF  observations  at  Millstone  Hill  showed  that  angle  scintil- 
lation was  generally  much  more  pronounced  in  traverse  angle  than  in  elevation 
angle.  This  is  presumably  related  to  the  high  inclination  angle  of  the  earth's 
magnetic  field  and  the  fact  that  the  irregularities  are  elongated  in  the 
direction  of  the  field. 
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The  rms  fluctuation  in  traverse  angle  about  the  mean  over  a 

three-second,  interval  (a  ) was  used  to  measure  angle  scintillation  and  the 

IK 

Millstone  Hill  study  served  to  establish  some  relationships  between  the  a 

IK 

and  the  amplitude  scintillation  index  S at  UHF.  Figure  3^  shows  a v S 

IK 

during  a period  of  intense  scintillation  occurring  in  a magnetic  storm.  It 

can  be  seen  that  the  mean  rms  angular  scintillation  increases  with  S,  but 

the  two  are  not  uniquely  related.  Whereas  S is  necessarily  limited  to  a 

saturation  value  near  unity,  the  angle  fluctuations  apparently  continue  to 

increase  as  the  irregularities  grow  more  intense. 

The  relationship  between  ct™  and  S using  all  satellite  tracks 

IK 

is  shown  in  Figure  35-  In  this  analysis,  the  scintillation  index  was  divided 

into  a number  of  ranges  and  for  each  range  the  distribution  of  a values 

IK 

was  constructed.  The  mean  and  most  probable  value  for  each  distribution 

were  calculated  and  are  given  in  Figure  35- 

The  scintillation  model  for  tracking  simulation  studies  described 

above  may  be  extended  to  cover  traverse  angle  fluctuations  by  using  the 

model  value  of  S obtained  from  equation  22  to  predict  a model  value  of  o™ 

IK 

from  Figure  35*  Before  this  procedure  is  satisfactory  some  thought  must  be 
given  as  to  how  the  results  in  Figure  35  may  be  generalized  to  all  frequencies 
rather  than  just  the  observing  frequency  of  400  MHz. 

The  frequency  scaling  law  for  angle  scintillation  needs  to  be 
established  before  any  of  the  Millstone  Hill  observations  made  at  400  MHz  can 

be  used  to  estimate  effects  at  higher  frequencies.  A detailed  analysis  of 

7 

the  spectra  of  the  amplitude  fluctuations  observed  at  Millstone  has  shown 
that  the  irregularities  in  the  ionosphere  that  are  responsible  have  a power - 
law  scale-size  distribution  with  an  index  of  four  for  the  three-dimensional 
wave  number  spectrum.  Theoretical  analysis  of  the  propagation  of  radio  waves 
through  a physically  thick  screen  containing  such  irregularities  has  shown 

that  for  weak  scattering  the  scintillation  index  (s)  will  scale  with  frequency 

_2  /q  -2 

( f ) as  f 7 , while  the  rms  angle  fluctuations  scale  as  f . Evidently, 

at  higher  frequencies  the  angle  fluctuations  will  rapidly  become  insignificant 

even  when  amplitude  fluctuations  can  be  detected. 
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SCINTILLATION  INDEX  S 

Fig.  34  RMS  angular  wander  of  target  position  in  traverse  at  UKF  vs  the 
amplitude  scintillation  index  S observed  during  a magnetic  storm. 
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Fig.  35  The  mean  and  most  probable  values  of  as  a function  of  scintillation 
index  from  UHF  satellite  beacon  tracking  records  gathered  at  Millstone  Hill  from 
January  1971  to  March  1973 • 
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As  there  is  an  approximate  linear  relationship  between  and. 

S for  S < .5  the  frequency  scaling  laws  suggest  that  this  portion  of  the 

Figure  35  can  be  generalized  for  any  frequency  ( f in  MHz)  of  observation  by 

adopting  a™  (-vv-)  as  the  ordinate  rather  than  <j  . For  larger  values  of  S 
IK  I IK 

saturation  effects  are  evident  and  in  any  case  the  frequency  scaling  laws 
become  uncertain  because  the  quoted  relationships  are  only  applicable  for 
weak  scatter.  In  view  of  these  uncertainties  the  best  approach  is  probably 
to  assume  that  the  figure  as  a whole  can  be  applied  at  any  frequency  by 
using  (“jr1)  as  an  ordinate  rather  than  . 

The  probability  distribution  of  traverse  angle  fluctuations  may 
be  taken  as  Gaussian,  while  the  correlation  time  may  be  assumed  to  be  the 
same  as  for  amplitude  fluctuations  (equation  30).  No  adjustment  to  the 
traverse  angle  fluctuation  is  necessary  to  allow  for  the  two-way  path 
involved  in  a radar  situation. 

K.  Model  for  Severe  Ionospheric  Scintillation 

A crude  model  of  severe  scintillation  effects  has  been  constructed 
to  specify  the  upper  limit  to  the  magnitude  of  the  fluctuations  that  might 
be  expected  for  a radar  operating  at  L-band  (1000  MHz).  The  maximum  degrada- 
tion of  radar  performance  caused  by  ionospheric  scintillation  can  be  deter- 
mined by  using  this  model  in  tracking  simulation  exercises.  During  geomag - 
netically  disturbed  conditions  the  probability  of  occurrence  and  intensity 
of  scintillation  increases  markedly  and  a Hrarse-case  model'1  was  based  on 
UHF  measurements  made  at  Millstone  Hill  at  a period  of  severe  magnetic 
disturbance  (K  = 8+).  During  this  time  the  rms  traverse  angle  fluctuation 
reached  the  largest  value  observed  during  the  2y  years  of  the  Millstone 
satellite  tracking  program , namely , 120  mdeg  (Figure  3^).  The  scintillation 
index  at  UHF  reached  the  saturation  level  near  unity  and  thus  could  not  be 
used  directly  as  a basis  for  extrapolation  to  L-band.  The  traverse  angle 

rms  did  not  exhibit  any  apparent  saturation  effects  and  scaling  was  accom- 

_2 

plished  by  assuming  the  f scaling  law  for  angle  fluctuations,  giving  an 

estimate  of  c % 20  mdeg  at  L-band  (1000  MHz).  Although  this  scaling  law 
TK 

is  supported  both  theoretically  and  experimentally  under  conditions  of  weak 
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scintillation,  its  use  in  the  present  case  must  be  considered  as  a reasonable 
guess  only.  Using  the  curve  in  Figure  35  which  relates  the  mean  value  of 

hoo 

f 

of  S = 0.55  for  the  worse-case  model  at  L-band. 

L.  Summary  of  Scintillation  Model 

1.  Amplitude  Scintillation 

A mean  scintillation  index  S (A)  as  a function  of  invariant 
latitude  is  obtained  from  the  appropriate  one  of  Figures  30a  to  30f.  The 
scintillation  index  for  the  observing  frequency  and  target  location  is  then 
estimated  from  equation  22  and  converted  to  a rms  amplitude  fluctuation  a 

R 

using  equation  26.  For  the  worse-case  model  at  f=  1000  MHz  a model  value 

of  <j  = .27  is  adopted  independent  of  target  location , provided  it  is  above 
R 

the  irregularities.  Based  on  a > a value  of  the  fluctuation  in  relative 

R 

amplitude  R for  a one-way  transmission  is  selected  from  a distribution  of 

_2 

the  form  given  by  equation  27>  with  m = (2a).  The  actual  radar  amplitude 
is  then  obtained  using  equation  29*  The  correlation  time  may  be  calculated 
from  equations  30  and  31>  using  the  values  of  E and  T for  the  track.  The 
radar  amplitude  should  be  held  constant  for  a time  t before  selecting  a new 
value  from  the  model. 

2 . Angle  Scintillation 

The  rms  fluctuation  in  traverse  angle  calculated  from 

the  model  value  of  S through  Figure  35*  For  the  worse -case  model  at  f =1000 

MHz,  <7^  is  taken  as  20  mdeg.  The  actual  traverse  angle  fluctuation  is  selec- 
1R 

ted  from  a Gaussian  distribution  with  zero  mean  and  standard  deviation  a . 

iR 

This  value  should  be  held  for  a time  t,  as  for  the  amplitude. 

M.  Depolar i zation 

A powerful  discriminate  against  false  targets  is  based  on  the 
depolarization  of  circularly  polarized  signals.  There  is  concern  that, 
during  disturbed  conditions,  the  auroral  ionosphere  could  introduce  suffi- 
cient depolarization  to  render  this  discriminate  useless.  This  possibility 
was  tested  by  R.  Crane  by  analyzing  the  UHF  beacon  tracking  results  obtained 
in  the  Millstone  Hill  propagation  study  (Section  II ). 


* = 


= 13  mdeg,  provides  an  estimate 


TR 


to  S,  and  an  ordinate  value  of  a 


TR* 
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In  many  cases  these  signals  exhibited  severe  angular  and 
amplitude  scintillation  along  paths  at  elevations  below  about  20-30°  to  the 
north  of  Millstone.  Simultaneous  observations  of  both  left-  and  right-hand 
circular  polarizations  were  recorded.  (The  transmissions  were  nominally 
right-hand  circular,  but  proved  in  practice  to  be  elliptic ally  polarized.) 

The  polarization  state  changed  slowly  with  changes  in  satellite/receiver- 
station  geometry.  The  gain  of  the  orthogonal  polarization  receiver  was 
controlled  by  the  primary  polarization  AGC  signal.  The  AGC  control  system 
was  effective  in  removing  fluctuations  of  limited  dynamic  range  that  occurred 
simultaneously  on  both  channels  at  frequencies  up  to  250  Hz.  For  strong 
fluctuations  with  peak-to-peak  spreads  of  more  than  20  cLB,  the  AGC  system 
did  not  remove  all  of  the  simultaneous  fluctuations  from  the  orthogonal 
channel  output,  and  residual  fluctuations  could  be  detected.  For  weak 
scintillation,  only  fluctuations  on  the  orthogonal  channel  which  were  not 
correlated  with  the  principal  polarization  fluctuations  would  be  detected. 

The  August  1972  records  for  the  most  disturbed  tracks  observed 
during  the  three-year  period,  March  1970  - March  1973 > were  analyzed  to 
measure  the  correlation  between  the  right-  and  left-hand  circular  components. 
Figure  36  shows  the  rms  log  of  the  received  power  at  150  MHz  and  400  MHz  for 
the  principal  channel  during  a pass  near  0400  hours  GMT  on  4 August  1972.  It 
can  be  seen  that  at  six  minutes  the  UHF  signal  scintillation  became  fully 
developed. 

The  rms  variation  of  the  log  of  the  orthogonal  cnannel  amplitude 
and  the  correlation  coefficient  between  the  log  of  the  orthogonal  channel 
output  and  the  log  of  the  principal  channel  output  is  shown  in  Figure  37 
for  the  entire  pass.  For  weak  scintillation,  (i.e.,  before  six  minutes)  the 
output  is  near  receiver  noise  and  no  correlation  is  evident.  After  seven 
minutes,  the  scintillation  is  much  stronger  and  a low  level  variation  is 
evident  in  the  orthogonal  channel  data.  However,  this  residual  output  is 
highly  correlated  with  the  scintillation  in  the  principal  polarization 
channel.  The  data,  therefore,  show  no  uncorrelated  fluctuations. 
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Fig.  36  RMS  log  (signal)  vs  time  at  150  and  400  MHz  for  a very  disturbed 
satellite  pass  in  August  1972.  It  can  be  seen  that  the  150  MHz  signal  amplitude 
scintillation  becomes  completely  developed  before  that  at  UHF . 
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Fig.  37  RMS  log  (signal)  vs  time  of  the  orthogonally  polarized  signal  com- 
ponent and  the  correlation  coefficient  between  the  fluctuations  at  VHF  and  UHF, 
for  the  pass  shown  in  Figure  36.  At  low  scintillation  and  AGC  circuit  removes 
the  fluctuations,  but  fails  to  do  so  for  strong  scintillation.  Fades  at  both 
polarizations  are  found  to  be  highly  correlated. 
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These  observations  show  that  the  fluctuations  for  both  polari- 
zations are  correlated,  i.e.,  both  polarizations  fade  together.  Thus,  though 
the  fading  itself  may  be  a nuisance,  discrimination  algorithms  that  depend 
upon  observations  of  the  depolarization  of  circularly  polarized  waves  will  not 
be  affected  adversely. 

VIII.  SUMMARY 

Table  1 attempts  to  summarize  the  maximum  values  of  the  radar 
metric  errors  discussed  above  that  are  liable  to  be  encountered  at  1000  MHz 
at  low  elevations  at  a location  near  the  auroral  zone.  It  should  be  under- 
stood that  the  errors  will  generally  be  less  than  given  here  as  the  reported 
values  in  most  cases  were  obtained  from  the  worse  cases  observed  over  the 
duration  of  the  Millstone  Hill  propagation  study. 

In  particular,  the  scintillation  result  represents  conditions  during 
a severe  magnetic  storm  and  there  must  be  some  doubt  as  to  the  appropriate 
frequency  scaling  laws  when  such  strong  scintillation  was  present  at  UHF. 
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TABLE  I 


SOURCE  OF  ERROR 


FREQUENCY 

MAXIMUM  ERROR  DEPENDENCE 


Effects  in  the  Ambient  Ionosphere 
Elevation  Bias 
Range  Bias 

Trough  Effects 

Elevation  Error 
Range  Error 


+ 7 radeg 
+ 75  m 

+ 9 mdeg 
- 10  m 


TIPs 

X>  100  km 

Elevation  Error  (peak-to-peak)  10  mdeg 
Range  Error  (peak-to-peak)  5 m 
X<  100  km 

Elevation  Error  (peak-to-peak)  10  mdeg 
Scintillation 

Relative  amplitude  (rms)  0.3 

Azimuth  (rms)  20  mdeg 


f 

f 


-2 

-2 


f 

f 


-2 

-2 


f 

f 


-2 

-2 


f 


-2 


f 

f 


-3/2 

-2 
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